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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  develop  a  coupled  ocean-estuary  simulation  tool  that  will 
consist  of  the  Regional  Oceanic  Modeling  System  (ROMS)  and  the  nonhydrostatic,  unstructured-grid 
SUNTANS  model,  in  order  to  study  the  interaction  of  coastal  regions  with  estuaries.  The  scientific 
goal  will  be  to  understand  the  dynamics  of  ocean-estuary  exchange  between  San  Francisco  Bay  and  the 
coastal  Pacific  Ocean,  with  ROMS  simulating  the  coastal  Pacific  Ocean,  and  SUNTANS  simulating 
San  Francisco  Bay. 

OBJECTIVES 

The  primary  objective  of  this  project  is  to  develop  and  implement  a  coupled  hydrodynamics  and 
sediment-transport  simulation  tool  to  study  the  complex  three-dimensional,  coupled  ocean-estuary 
processes  that  influence  exchange  through  the  Golden  Gate  channel  in  San  Francisco  Bay.  In  order  to 
understand  the  influence  of  regional  processes  in  the  coastal  ocean  and  how  they  are  coupled  to  local 
processes  in  San  Francisco  Bay,  the  project  will  employ  a  one-way  coupled  simulation  tool,  in  which  a 
model  for  San  Francisco  Bay  is  coupled  to  a  regional  model  for  the  coastal  Pacific  Ocean.  The 
SUNTANS  model  (Fringer  et  ah,  2006a)  will  be  employed  to  simulate  flow  in  San  Francisco  Bay,  and 
it  will  obtain  boundary  conditions  at  the  coastal  ocean  from  simulations  of  the  California  coastal 
current  using  ROMS  (Shchepetkin  and  McWilliams,  2005).  SUNTANS  is  well-suited  to  simulate  the 
exchange  flow  at  the  Golden  Gate  because  it  employs  an  unstructured  grid,  which  enables  simulation 
domains  that  include  the  entire  Bay  while  also  resolving  finescale  flow  features  near  the  Golden  Gate. 
Furthermore,  because  SUNTANS  is  a  nonhydrostatic  code,  it  will  capture  nonhydrostatic  features  of 
the  exchange  flow  which  may  be  important  at  high  resolution.  ROMS,  on  the  other  hand,  was 
originally  written  to  simulate  regional  flows  like  the  California  Coastal  Current,  and  employs  highly 
accurate  methods  to  simulate  the  complex  dynamics  of  upwelling  fronts,  jets,  and  squirts  over  seasonal 
time  scales  (Marchesiello,  et  ah,  2003).  Using  a  stand-alone  SUNTANS,  simulations  will  be 
performed  in  idealized  domains  to  understand  the  fundamental  mechanisms  governing  three- 
dimensional  exchange.  In  the  real,  field-scale  domain,  one-way  coupling  will  be  used  to  understand 
how  tidal  and  seasonal  variability  in  the  coastal  currents  as  well  as  in  the  Bay  influence  the  dynamics 
of  ocean-estuary  exchange  through  the  Golden  Gate  channel. 
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APPROACH 


Simulations  of  San  Francisco  Bay  will  employ  SUNTANS  (Stanford  Unstructured  Nonhydrostatic 
Terrain- following  Adaptive  Navier-Stokes  Simulator),  a  free-surface,  nonhydro  static,  unstructured- 
grid,  parallel  coastal  ocean  and  estuary  simulation  tool  that  solves  the  Navier-Stokes  equations  under 
the  Boussinesq  approximation  (Fringer  et  al., 2006a).  The  formulation  is  based  on  the  method  outlined 
by  Casulli  &  Walters  (2000),  in  which  the  free-surface  and  vertical  diffusion  are  discretized  with  the 
theta-method,  which  eliminates  the  Courant  condition  associated  with  fast  free-surface  waves  and  the 
elevated  friction  term  associated  with  small  vertical  grid  spacing  at  the  free-surface  and  bottom 
boundary.  For  flows  with  extensive  wetting  and  drying  such  as  San  Francisco  Bay,  advection  of 
momentum  is  accomplished  with  the  semi-Lagrangian  advection  scheme  (Casulli  and  Cheng,  1992; 
Staniforth  and  Temperton,  1986),  which  ensures  stability  in  the  presence  of  cells  that  fill  and  empty 
with  the  tides. 

To  simulate  the  influence  of  coastal  currents  on  exchange  between  San  Francisco  Bay  and  the  coastal 
Pacific,  the  SUNTANS  simulations  will  be  nested  within  the  ROMS  model  (Shchepetkin  and 
McWilliams,  2005),  which  also  has  a  nonhydrostatic  module  (Kanarska,  et  al.,  2007).  ROMS  is 
discretized  in  horizontal  curvilinear  coordinates  and  has  a  generalized  terrain-following  vertical 
coordinate  that  is  configured  to  enhance  resolution  near  the  sea  surface,  and  the  time-stepping 
algorithm  is  split-explicit. 

The  simulations  will  employ  a  one-way  coupled  ROMS-SUNTANS  tool  that  has  been  developed  and 
applied  to  simulate  nonlinear  internal  tides  on  the  Australian  North  West  Shelf  (Fringer  et  al.,  2006b), 
and  development  of  the  associated  two-way  coupled  tool  is  underway  in  collaboration  with  Prof.  J. 
McWilliams  at  UCLA  in  order  to  study  tide-current  interactions  in  Monterey  Bay  as  well  as  the  South 
China  Sea.  The  SUNTANS  simulations  will  be  one-way  nested  within  simulations  of  the  California 
Coastal  Current  (Marchesiello  et  al.,  2003).  These  simulations  will  provide  sea-surface  height, 
velocity  fields,  and  sediment,  salinity  and  temperature  at  the  offshore  boundary  for  SUNTANS. 
Because  the  ROMS  simulations  of  the  coastal  current  are  nontidal,  tidal  currents  from  the  OTIS  tidal 
model  (Egbert  et  al.,  1994)  are  included  by  adding  them  to  the  mesoscale  currents  provided  by  ROMS 
at  the  SUNTANS  boundaries.  In  one-way  nesting,  the  boundary  conditions  supplied  by  ROMS 
influence  SUNTANS,  yet  it  is  assumed  that  the  flowfield  as  computed  by  SUNTANS  does  not 
influence  the  ROMS  simulation  (see,  e.g.  Fringer  et  al.,  2006c).  The  boundaries  of  the  SUNTANS 
simulations  will  be  placed  in  order  to  minimize  their  potential  influence  on  the  exchange  plume  as  it 
evolves  on  the  shelf  offshore  of  the  Golden  Gate.  As  long  as  the  plume  that  exits  the  SUNTANS 
boundary  does  not  have  the  potential  to  reenter  the  Bay,  then  the  boundary  is  far  enough  from  the 
mouth  of  the  Bay  and  one-way  nesting  is  justified. 

Because  both  SUNTANS  and  ROMS  are  well-established  models  and  the  one-way  coupling 
framework  has  already  been  developed  (Fringer  et  al.,  2006b),  a  bulk  of  this  project  will  focus  on 
implementation  of  the  models  to  understand  detailed  physics  of  ocean-estuary  exchange.  Yet  while 
SUNTANS  has  successfully  been  applied  to  study  an  estuary  with  significant  wetting  and  drying 
(Wang  et  al.,  2008),  advection  of  momentum  in  the  presence  of  wetting  and  drying  presents  significant 
challenges  when  high  grid  resolution  is  employed.  While  accuracy  of  semi-Lagrangian  advection 
schemes  depends  on  the  time-accuracy  of  the  traceback  as  well  as  the  spatial  accuracy  of  the 
interpolation  (Staniforth  and  Cote,  1991),  typical  implementations  of  semi-Lagrangian  advection  for 
estuarine  simulations  with  coarse  resolution  employ  first-order  accurate  spatial  interpolation  (Casulli 
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and  Cattani,  1994)  and  ignore  vertical  advection  of  momentum  (Zhao,  2007).  While  this  does  not 
incur  problems  when  coarse  resolution  is  employed,  high-resolution  requires  implementation  of 
vertical  advection  and  higher-order  interpolation  in  the  presence  of  complex  bathymetry.  This  is  not  a 
trivial  task,  since  it  is  difficult  to  ensure  stability  in  the  implementation  of  higher-order  interpolation 
schemes  on  unstructured,  finite-volume  grids,  particularly  with  wetting  and  drying  (Zhao,  2007). 
However,  if  the  stencil  is  suitably  chosen,  stability  can  be  ensured,  as  pointed  out  by  Leonard  (2002). 
In  the  proposed  work  we  will  explore  stencils  for  the  Lagrangian  traceback  that  improve  stability  when 
second-order  accurate  interpolation  is  used. 

WORK  COMPLETED 

Funding  began  mid  FY08  so  no  work  has  been  completed  to  date. 

IMPACT/APPLICATIONS 

The  coupled  simulation  tool  that  will  be  developed  for  this  project  has  the  potential  to  form  the  basis 
for  a  wide  range  of  modeling  studies  of  other  important  coupled  ocean-estuary  systems. 

RELATED  PROJECTS 

The  California  Coastal  Conservancy  has  recently  funded  the  development  of  an  open-source 
hydrodynamics  and  sediment  transport  model  using  SUNTANS  to  assess  the  impacts  of  marsh 
restoration  projects  on  salt  and  sediment  dynamics  in  San  Francisco  Bay.  Fringer  is  involved  in  this 
project  in  collaboration  with  Profs.  Mark  Stacey  and  Zack  Powell  at  U.  C.  Berkeley  and  Profs.  Jeffrey 
Koseff  and  Stephen  Monismith  at  Stanford  University. 
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